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We study the ground state energy Eg{u) of N classical 0(n) vector spins with the Hamiltonian 
TL — — Si>j JijSi.Sj where the coupling constants {Jij} are arbitrary. We prove that So(n) is 
independent of n for all n > n ma x{N) = I vgjSI^ i j s ] low that this bound is the best possible. 
We also derive an upper bound for Ea(m) in terms of Ea(n), for m < n. We obtain an upper bound 
on the frustration in the system, as measured by F(n) = ^*>jJ J ' 3 _l££lli _ "We describe a procedure 
for constructing a set of Jij's such that an arbitrary given state, {Si}, is the ground state. 
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I. INTRODUCTION 

In this paper, we study the ground states of N unit 
classical 0(n) spins, Si, having a hamiltonian of the form 



H = 



JijSi.Sj 

i>j 



(1) 



where Jij's are arbitrary real numbers-positive, negative 
or zero. Such hamiltonians with arbitrary bonds and cou- 
plings are of interest in the context of disordered systems, 
especially spin glasses [1]. One of the interesting ques- 
tions is the behavior of the ground state energy as the 
spin space dimension, n, is increased. For example one 
can study the behavior of such models when n is large. In 
this context, Hastings [2] proved that for N spins beyond 
a spin space dimension of n max {N) = [ V 8jv + 1 ~ lj the 
ground state energy does not decrease any further and 
also that this bound is saturated. Aspelmeier and Moore 
[3] have then used this bound in accelerating their numer- 
ical simulations of spin glasses. We provide an alternate 
proof for this bound. A similar analysis has been done 
earlier in the context of correlation matrices by Grone, 
et. al. [4] and for a relaxed version of maxcut problem 
of theoretical computer science [5]. 

An interesting question is the behavior of the average 
number of non-zero spin space components (average over 
disorder) in the ground state, as a function of the num- 
ber of spins, N. For the infinite range model with gaus- 
sian distributed Jy's, this number increases as where 
fi — | [2,3]. Lee, Dhar and Young [6] have numerically 
determined fj, for several different models. 

We also derive both upper and lower bounds on the 
ground state energy of 0(m) spins in terms of the ground 
state energy when they are replaced by 0(n) spins (to ^ 
n) keeping the couplings, Jij's, the same. A stronger 
bound is also provided for Ising spins (m = 1) when all 
couplings are antiferromagnetic and Ec(n) is low. 

We also consider the problem of finding the ground 
state of such a hamiltonian [7]. We study the inverse 



problem — how to find a (non-trivial) hamiltonian of the 
form in equation (1) so that a given spin state {Si} is the 
ground state. This question is trivial for Ising spins. One 
just assigns a non-negative Jij if the spins are parallel, 
and a non-positive Jij if they are anti-parallel. However 
with 0(n) spins (n > 1) the problem is non-trivial and in 
some cases there is no (non-zero) solution, for example, 
A = 3, n > 3, with the three spins non-coplanar. In gen- 
eral, to find the desired set of couplings, Jij's, we can ex- 
press the hamiltonian in terms of the angles of the spher- 
ical polar coordinates of the spins and set the derivatives 
with respect to the angles equal to zero at the angles cor- 
responding to the desired ground state. This gives a set of 
linear relations between the couplings, Jij's. In addition, 
to ensure that this extremum is a minimum, and not a 
maximum or a saddle point, we have the additional con- 
straint of the Hessian being positive semi-definite. Find- 
ing Jij 's which simultaneously satisfy the linear relations 
as well as the positive semidefiniteness constraint on the 
Hessian is non-trivial. It is a semi-definite programming 
problem [8] for which fast algorithms and their software 
implementations are available. 

We provide a simple procedure for obtaining a large 
class of such hamiltonians. However, not all hamiltoni- 
ans with {Si} as the ground state are obtained by this 
procedure. We conjecture a characterisation of the hamil- 
tonians obtained and give proof of a part of the conjec- 
ture. 

The plan of the paper is as follows. In section II 
we summarise some properties of the correlation ma- 
trices of classical spin states which are found useful in 
the later sections. In section III, we prove that for N 
spins beyond a spin space dimension of n max (N) the 
ground state energy becomes independent of n. That 
this is the best bound is proved by providing a sequence 
of graphs and couplings, one for each N, such that 
Eaijimax - 1) > EG{n max ). In section IV, we derive 
both upper and lower bounds on the ground state energy 
of 0(m) spins in terms of the ground state energy when 
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they are replaced by 0(n) spins (m ^ n) keeping the cou- 
plings, Ji/s, the same. For Ising spins with all couplings 
antifcrromagnetic, in a special case, a stronger bound is 
derived. We obtain an upper bound on spin frustration, 
as measured by F(n) = ^'>jJ l7 ' J l+- E j 3 ("- ) — we show that 

F(n) — F(oo) < (3 n where /3„ is a constant, independent 
of Jy's. In section V, we provide a procedure for con- 
structing hamiltonians of the form in equation (1) with 
arbitrary given state {Si} as the ground state. Section 
VI summarises the results. 



II. SOME PROPERTIES OF CORRELATION 
MATRICES OF CLASSICAL SPIN STATES 

For an arbitrary state {Si} of 0(n) spins define the 
(N x N) correlation matrix, C = [Si.Sj]. Alternatively, 



C^S'S 



(2) 



th 



where S is the (n x N) matrix with vector of the i 
spin as the i th column. Clearly, C is real, symmetric, 
has diagonal elements unity and can be written as C = 
ODO T , where O is an orthogonal matrix and D diagonal. 

C is positive semidefinite, i.e. all eigenvalues of C 
are non-negative, since for every x £ R N ,x T Cx = 
(Sx) T (Sx) > 0. 

The number of non-zero (and hence positive) eigenval- 
ues of C is atmost n. This can be seen as follows: each 
row of C is a linear combination of the n rows of S im- 
plying that the number of linearly independent rows of C 
is atmost n. Diagonalising C, let C — ODO T where O is 
an orthogonal matrix and D diagonal with (let's say) first 
k eigenvalues positive and rest zero. The rows (columns) 
of O arc mutually orthogonal and hence linearly inde- 
pendent. (DO T ) now has k linearly independent rows 
and thus ODO T also has k linearly independent rows. 
The number of linearly independent rows of C we have 
already argued to be atmost n. Hence the number of 
positive eigenvalues of C is atmost n. 

Conversely, if C is a real, symmetric matrix with di- 
agonal elements unity and having n or fewer positive 
eigenvalues and rest zero, then there exists a spin state 
of classical 0{n) unit spins for which it is the correla- 
tion matrix. To see this, C being real and symmetric, 
can be diagonalised as C = ODO T = {OVD){OVD) t 
where O is an orthogonal matrix and D diagonal with 
first k(< n) diagonal entries positive. The last (N — n) 
rows of (0\fD) T are known to be zero and we drop them 
to define an (n x N) matrix S such that C = S T S. Since 
cu = 1 V i each column of S can be interpreted as a unit 
classical 0(n) spin, C being their correlation matrix. 



III. THE INDEPENDENCE OF THE GROUND 
STATE ENERGY FROM n FOR n > n max {N) 

Consider the variation of the ground state energy 
Ec(n) as a function of the spin space dimension n of the 
0(n) spins keeping the couplings, J^'s, the same. For 
n' > n we have Ea(n) > Ea(n') because for any state 
of 0(n) spins we can construct a corresponding state of 
O(n') spins with the same value of energy by augmenting 
each vector with (n' — n) zeroes. Also for any n > N we 
have Ec{n) = Eq{N) because N spins span an atmost 
TV dimensional subspace of the n dimensional spin space 
implying that by an appropriate choice of basis we can 
make all coordinates after the first N coordinates zero 
and by dropping them we get an O(N) spin state with 
the same value of energy. 

Theorem 1 For N classical unit 0(n) vector spins 
with hamiltonian H — — X^>j JijSi-Sj, where {Jij} are 
any real numbers, the ground state energy Eq(ji) = 
E G (n m ax) for all n > n max where 



m = [ 



VSN + 1 - 1 



(3) 



Here [x\ for x e R is the greatest integer not greater 
than x. 

Proof Let us summarise the idea of the proof be- 
fore getting into the details. Suppose we have a ground 
state which has more than n max (N) dimensions. Start- 
ing fron the correlation matrix of this state we con- 
struct another matrix which is the correlation matrix of 
a spin state which is embedded in one less spin space 
dimension but has the same energy. This construc- 
tion always works whenever the spin state is embed- 
ded in more than n max {N) dimensions. Since, as shown 
above, Ec(n) > Ea(n') for n' > n this implies that for 
n > n max , E G (n) = E G (n max ), as desired. 

We now discuss the proof. When the spins are O(N) 
vectors let {Si} be a ground state. Consider the cor- 
relation matrix, C, with elements dj — Si.Sj for all i, 
j. Diagonalising C, we can write C — T DO with D a 
diagonal matrix and O an orthogonal matrix. Let 



D = 



dl . 


. 


.. 


. 






.. 


. 


. 


■ d k 


.. 


. 


. 


. 








. 










(4) 



where the first k diagonal entries of D are positive and 
rest zero. 

Consider C = O t (D + rB)0 where B is symmetric 
with Bij = if i > k or j > k. This leaves \k{k + 1) free 
parameters in B and ensures that C is also symmet- 
ric and the zero eigenvalues of C and the corresponding 
eigenvectors are not perturbed. Also, let B satisfy 
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[0 T B0] tt = 



(5) 



for all t = 1,2, ...N. This ensures that the diagonal 
elements of C remain unchanged. 

The ^k(k + 1) free parameters of B must satisfy the 
N linear homogenous equations (5). Hence whenever 
^k(k + 1) > N such a non-zero B will exist and we can 
increase r till one of the first k eigenvalues of C becomes 
zero. Thus we obtain a matrix C which is the correlation 
matrix of a spin state embedded in (k— 1) dimensions. As 
shown in the next paragraph, this spin state is a ground 
state. Thus applying this procedure repeatedly we obtain 
a ground state embedded in atmost [ ^ /8jv ^ 1 ~ 1 j dimen- 
sions. 

The matrix B chosen above is such that for r small 
enough C = C ± rO T BO are both correlation ma- 
trices of valid spin states with energy — J2i>j Jij c ij ^ 
r Ei>j . Jij(O r BO)ij]. Since we started from a C 
which was a ground state this can happen only if 
T,i>j ■h ] {0 T BO) l3 = i.e. if T BO was a neutral di- 
rection. Thus the correlation matrix C also corresponds 
to a ground state. 

Hence we have provided a construction for continu- 
ously deforming a ground state and bringing it to lie in an 
atmost n max (N) dimensional subspace of the spin space 
without changing the energy, thus proving the desired 
result. 

Theorem 2 The bound in theorem 1 is the best possi- 
ble, i.e. there exist values of {Jij} such that Ec{n max ) < 
E G {n max - 1) where n max (N) = [ v^v+T-i j . 

Proof Consider three spins S p , S q and S pq with J pq = 
—J and Jp(p q ) — Jq( P q) — V%J, see figure 1. It is easy 
to see that for this system of three spins in the ground 
state S p is perpendicular to S q . If wc integrate over S pq , 
we get an effective interaction between S p and S q . Now 
construct a set of k spins with this effective interaction 
between every pair of them (including the intermediate 
spins the total number of spins will be N — ^k(k + 1)). 
Their energy gets fully minimised only when these k spins 
are perpendicular to each other which happens only in an 
atleast n max = k dimensional space thus completing the 
proof. 



IV. BOUNDS ON THE GROUND STATE 
ENERGY 

We have seen that Ec(m) > Ec(n) for m < n. Now 
we will derive an upper bound on Ecim) in terms of 
Ec{n). This result (theorem 3) generalises a known re- 
sult on the performance of Goemans- Williamson algo- 
rithm for maxcut problem of theoretical computer sci- 
ence [5]. The result by Goemans and Williamson, when 



J4=6 N=10 N=(l/2)k(k+l) 

S ^fT^\ s Complete graph 
of k spins where 
each edge is itself 
an S subgraph 

5 2 



An S 2 subgraph between two spins Sp and means: 





FIG. 1: A sequence of examples for which Ea{n ma x(N)) < 
Ec(n m ax{N) — 1). Shown are three of those members of 
the sequence for which V8W+1-1 j s an integer (such AT's are 
3, 6, 10, ... = ifc(fc+l)). Rest of the members are obtained by 
adding appropriate number of free spins to the example of the 
last N for which v^siv+i-i ; s an integer eg. N = 8 example 
has two more free spins added to the N = 6 example. 



translated into statistical physics language, would cor- 
respond to the special case of m = 1. Theorem 4 is a 
translation of a known result on maxcut problem into 
statistical physics language [5]. The connection between 
the problem of finding the ground states of Ising spins 
and maxcut problem has been known before [9]. 

It will be helpful to summarise the overall strategy 
before getting into the details. Suppose the various pos- 
sible orientations of 0(m) spins occur according to an 
arbitrary given probability distribution. Then the en- 
ergy is also a random variable and the expected value of 
the energy will always be greater than or equal to the 
ground state energy, i.e. Ec(m) < E[H m ], where E[H m ] 
denotes the expected value of the energy of 0(m) spins. 
If we choose the probability distribution in such a way 
that we are able to bound E\H m \ in terms of Ecin) from 
above we would have obtained the desired result. 

Now we give the derivation in detail. First we define a 
randomised procedure for obtaining an 0(m) state, say 
{Si( m )}, from the ground state {S'^A of 0(n) spins. 
In the spin space of 0(n) spins randomly choose an Tri- 
dimensional subspace and project all the spins onto it. 
Normalise the O(m) vectors thus obtained. Clearly dif- 
ferent 0(m) states are obtained by this procedure de- 
pending on which m-dimensional subspace was chosen 
for projection. The expectation value of the 0(m) en- 
ergy is E[H m ] = - J2i>j Jij E \Si( m )-Sj( m )\- Now if V mn 
denotes a projection operator from n to m dimensions 



E[Si(m)-Sj(m)] = / — ^ TTT 5 -dVmn = fmn(0ij) 

where % is the angle between S^,-. and S'-, and the 
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integral is over all projection operators V mn with equal measure. As an example, in spherical polar coordinates, 



f 23(e) 



/•7T />27T 
J cf,2=0 J(t>l=( 



sin <j>2 (cos 8 — cos 0i sin 2 2 cos(0i — 0)) 
Att^/I — cos 2 0i sin 2 </> 2 ^/l — sin 2 </> 2 cos 2 (0i — 0) 



(7) 



By reversing the direction of 5^ n ^ we observe that 

/mn(7T — 0y) = — ,/mn(%) (8) 

Also ^Z™"^ > for all e (0,tt]. Hence we can find 
a lower bound on \^}p ' denoted by a mn , which gives 

/mn(^) < (1 - ami.) + <^mn COS 

* (9) 

Also replacing by (7r — 0) in this inequality we get 

- fmn{0) < (1 - "mn) - %« COS (10) 

For Jjj < 0, using (9), we get 

— Jij fmnifiij) 1— ~ (1 — Ctmn)Jij OtmnJij COS 

(ii) 

For J i:) > 0, using (10), we get 

Jij fmniPij) — (1 C*'mn)Jij ^mnJij COS 0ij (1^) 

Summing (11) over all those ij-pairs for which < 
and (12) over all those ij-pairs for which > and 
adding we get 

E\H m ] = - J~]Jijfmn{6ij) < (^-a mn )^\J i:j \+a mn E G (r 



S o 



Angle between the 0(4) spins (in radia 



FIG. 2: The functions / 34 (6») and 934(f) = \1^b 



i>j 



i>j 



(13) 

Now the minimum value of a random variable is always 
less than or equal to its expectaion value. 
Therefore, we have 
Theorem 3 For (m < n) 

E G {n) < E G {m) < (l-a ro „) ^ | Jy|+cwE G (n) (14) 

i>j 

where a mn is the minimum value of \^^p over the 
interval € (0, n] and / m „(0) has been defined above. 

Or rearranging the inequality, (^T - ) E G (m) — 

(^r)E i>j \Ji j \<E G (n)<E G (m) 
As an example, for m = 1 and n arbitrary, fi n (9ij) = 

■ 

1 — 2-^- and a\ n s=s 0.87856 [5j. We have determined / 2 3 
and /34 numerically by representing them as integrals in 
spherical polar coordinates, see equation (7) for instance 
The graphs of fa(0) and (734(0) 



As a specific instance, for the triangular lattice anti- 
ferromagnet E G (1) = —^J implying that — (0.52)J < 
E G {2) < -A J. It is known that E G (2) = -(0.5) J which 
compares very well with the non-trivial part of the in- 
equality. 

Now, we prove a stronger bound for a special case of 
Ising antiferromagnets: 

Theorem 4 For the special case of m = 1, let all 
couplings, Jij, be antifcrromagnetic. Then for the case 
(E G (n) < 8 J2i>j I) (<5 ~ —0.69) we have the stronger 
bound 



E G (n) 

< E G (1) 

< (-Ei^-i)l 

z — ; 7r 

i>3 



(15) 
(16) 



arccos 



E G (n) 



Proof Again using the randomised procedure in the 
derivation of theorem 3, from the ground state {<S^} 
of 0(n) model, various Ising states are obtained with 



different probabilities such that E[Hi] 



7T I 



where xa — S', v <S"v > (using 



Wgg) 

1— COS # 



are shown in 



figure 2. We find that a 23 ~ 0.96 and a 3 4 w 0.98. 



fln(0ij) = 1 

Consider the function arccos x. Draw the oblique tan 
gent from (1,0) to the curve, intersecting the curve tan- 
gentially at (8, arccos 8) . Consider the function h(x) 



■5 



which is the same as arccosa; for x < S and the same 
as the tangent for x G [<5, 1] . 

Clearly, arccosa;^ > h(xij) and since all J^ < 0, 



< 



— Jijh(xij) — Jij 



i>j 



i>j 



i>j \i>j 



P >q \Jpq\ 



(18) 
(19) 



where the last inequality uses the convexity of h(x). 

Since the minimum value of a random variable is less 
than or equal to its expectation value, we have Ec(n) < 
Eg{^) < E[Hi}. Also, for x < S we have h(x) — arccosa: 
and the desired inequality is proved. 

In the presence of antiferromagnetic Jij's, there may 
not exist any spin configuration that minimises the en- 
ergy of each individual bond to — | | . One of the possi- 
ble measures of the frustration of spins is 



F{n) = 



_ Ei^lJijl+Ecin) 



(21) 



We can consider spin frustration as arising in two steps: 
first we choose the Jij's but do not put any restriction 
on the dimensionality of the spin space — it is allowed 
to be as large as desired for the minimisation of energy. 
The frustration of this system will be F(oo) which will 
be the same as F(N) because Eo{n) = Eq(N) for n > 
N. To obtain the actual O(n) system we now restrict 
the number of dimensions in the spin space to n, thus 
increasing the spin frustration from F(oo) to F(n). 

Theorem 5 If (m < n) 



E G {m)~E G {n) 



< 2(1 



.) 



As a particular case, 

F(n)-F(oo) < 2(1 -a nN ) 



(22) 



(23) 



where a mn are the same as in theorem 3. 

Proof In theorem 3, subtract Ea{n) throughout, di- 



vide by J2i>j an d observe that ^ 



Eg{n) 



> -1 thus 



completing the proof. 

In particular, for all n, F{n) — F(oo) < F(l) — F(oo) < 
2(1 - auv) w 0.24288. 



V. PROCEDURE FOR CONSTRUCTING A 
MODEL WITH AN ARBITRARY GIVEN 
GROUND STATE 

For classical spins of 0(n) type let {S^} be a given 
state. We want to construct a hamiltonian with only 



two-spin Heisenberg type interactions which has {S^} as 
the ground state. The following procedure constructs a 
hamiltonian of the form in eq. (1) (upto a constant) 
which has {S^} as the ground state. 

1. For the given ground state {S'i} construct the corre- 
lation matrix, C , such that c^- = S^.Sj Vz, j = 1,2, ...N. 

2. Let C = 0'D'O lT where O' is an orthogonal matrix 
and D' diagonal with, let's say, the first k diagonal entries 
non-zero and rest entire matrix zero. 

3. Construct an (NxN) auxilliary matrix G as follows: 



G = 



G3 Gi 



(24) 



where G\ is a (k x k) matrix, etc. Moreover, choose 
Gi = 0, G 2 = 0, G 3 = and G 4 to be any (A - k) x 
(N — k) real, symmetric matrix with all eigen values non- 
positive. 

4. Define J = [Jy] = 0'GO' T . 

Theorem 6: For the hamiltonian Ti = 
~ Sfj=i JijSi-Sj thus constructed, the spin state 
{S'i} is the ground state. 

Proof: For any spin state {Si}, construct the correla- 
tion matrix G = [cy] = [Si-Sj] and diagonalise it, 



G = ODO T 



(25) 



where O is an orthogonal matrix and D is diagonal with 
all eigenvalues non-negative. 

Also for the matrix J defined above let 



(- j) T = 6bo T 



(26) 



where O is an orthogonal matrix and D is diagonal. Since 
J is negative semidefinite, (— J) is positive semidefinite, 
thus the entries of D arc non-negative. 

NowW = Tr((-J) T C). Using equations (25) and (26) 
and repeatedly using Tr(AB) = Tr(BA) we get H = 

Tr[W T W] > where W = \Zd6 t OVD- Therefore, for 
any state {Si}, 



H>0 



(27) 



im- 



For {Si}, by construction, H = Tr{(-J) T C) 
plying that {S^} is a ground state of H. 

Although a large number of hamiltonians with arbi- 
trary given state {S^} as the ground state can be ob- 
tained by this procedure, not all the hamiltonians with 
this property are obtained. For instance, it can be easily 
checked that for three Ising spins, one up and the other 
two down, happens to be a ground state when all three 
couplings are antiferromagnetic with equal strength, but 
this set of couplings can not be obtained by the above 
procedure for any allowed choice of the matrix G4. Thus 
we would like to characterise which hamiltonians can be 
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obtained by this procedure for a given ground state and 
which hamiltonians can not be obtained. 

We expect that a hamiltonian with {S^} as the ground 
state is obtained by this procedure if and only if upon 
replacing the given spins by spins with any higher spin 
space dimension, keeping Jy-'s the same, the ground state 
energy remains the same. The if part is our conjecture 
while the only if part is proved as follows: for any n > n, 
by augmenting each vector of {S'j} by (n — n) zeroes, 
we can obtain a state with the value of the hamiltonian 
Tin = 0. Since the hamiltonian is expressible as the trace 
of the product of two symmetric positive semidefinite ma- 
trices its value cannot be negative as in equation (27), im- 
plying that the n dimensional state thus obtained is the 
ground state of 0{n) spins. Therefore, for h > n we have 
Ec{ n ) — Eg(ji) thus completing the proof. This proof 
is consistent with the case of three Ising spins with an- 
tiferromagnetic couplings discussed above because if wc 
replace three Ising spins by XY-spins the ground state 
energy decreases from — | J per bond to — \ J per bond. 

VI. SUMMARY 

We showed that as we increase the spin space dimen- 
sion, n, the ground state energy, Ecin), becomes inde- 
pendent of n beyond a spin space dimension of n max (N) , 
and this bound is the best possible. For m < n we de- 
rived an upper bound for Eq^tti) in terms of Eq(u), the 
lower bound was trivial. A stronger version for a spe- 
cial case of rn = 1 was also proved. Similar bounds on 
Edm) — Ec{n) and a measure of spin frustration, F(n), 
were derived. A procedure was given for constructing a 
hamiltonian with an arbitrary given spin state, {S^}, as 
the ground state. 
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